
IETEJournal of Research 


ISSN: 0377-2063 (Print) 0974-780X (Online) Journal homepage: http://www.tandfonline.com/loi/tijr20 


Taylor &. Francis 

Taylor & Francis Croup 


Design of Roll Control System for Rockets by 
Simulation Techniques 

N. V. Kadam, S. Dasgupta, Kusum Chandragupta & N. Sundararajan 


To cite this article: N. V. Kadam, S. Dasgupta, Kusum Chandragupta & N. Sundararajan (1974) 
Design of Roll Control System for Rockets by Simulation Techniques, IETE Journal of Research, 
20:8, 381-385, DOI: 10.1080/03772063.1974.11487423 

To link to this article: http://dx.doi.org/10.1080/03772Q63.1974.11487423 



Published online: 21 Aug 2015. 



Submit your article to this journal GP 



View related articles GP 


Full Terms & Conditions of access and use can be found at 
http://www.tandfonline.com/action/journallnformation?journalCode=tijr20 


r 


Download by: [University of California, San Diego] 


Date: 24 June 2016, At: 23:21 


















Downloaded by [University of California, San Diego] at 23:21 24 June 2016 


Design of Roll Control System for Rockets by Simulation 

Techniques 1 

N. V. KADAM*, S. DASGUPTA : , KUSUM CHANDRAGUPTA* & N. SUNDARARAJAN* 

Manuscript received on 15 June 1973. 

This paper describes the application of simulation techniques both analog and digital 
to arrive at optimum control system parameters namely control thrust level, attitude and rate 
feedback gains for roll control of a sounding rocket. 


D URING an uncontrolled flight of a rocket 
through the atmosphere, an angular motion 
about its longitudinal axis (rolling motion) is 
developed due to aerodynamic disturbance torques 
arising out of fin-misalignment, thrust misalignment 
etc. In order to effect attitude control about the 
pitch and yaw axes of rockets during flight through 
the atmosphere, it becomes necessary to reduce this 
roiling motion to zero. The maintaining of roll attitude 
as small as possible calls for an active control system. 

This study describes the fixation of parameters 
for such a control system, by iheans of analog and 
digital simulation 1 ’ 2 . Though the design by simu¬ 
lation can be done for any rocket, it is illustrated 
here for a roll control system of a Rohini-125 sounding 
rocket. A preliminary estimation of parameters for 
roll control system for Rohini-300, assuming time 
invariant plant and constant disturbances has been 
given in (Tech, memorandum). 


SYSTEM DESCRIPTION 

The system is illustrated in the block diagram 
form in Fig. 1. The control law selected is of the 
‘Single level bang-bang control with dead zone* 
type. The activating function is obtained by using 
a combination of quantities proportional to roll angle 
and roll rate. 

The system dynamics about the roll axis is given by 
J t [/ (0 6(0 ] + MO = M) - u{t-t d ) (1) 

where I(t) — roll moment of inertia at time t 

6(0) 0(0 — the roll angle and roll rate res¬ 
pectively. 

A d (t) — Aerodynamic damping torque 
A d (t) = K d {t )0 where K d (t) is the aerodynamic 
damping coefficient. 

A(t )— Aerodynamic roll inducing torque 

u(t) — Control torque at t 

id — Control actuation delay time 


AERODYNAMIC 



Fig. 1. Block diagram 


tAn earlier version of this paper has been presented in March, 73 in a Symposium on Process Control, held at A.C.C.E.T., 
Karaikudi, Tamil Nadu, India. 

♦Control, Guidance and Instrumentation Division, Space Science and Technology Centre, Trivandrum. 





Downloaded by [University of California, San Diego] at 23:21 24 June 2016 


J. 1NSTN ELECTRONICS & TELECOM. ENGRS, Vol. 20, No 8,1974 


/(f) is assumed to be 

= [ 1(fo) ~ t /(4) ] 

for t < t b (2) 

= /(f 6 ) for t > t b 

where t 0 , 4 are the initial and bum out times respec¬ 
tively. The linear variation of moment of inertia 
is true only for end-burning type of propulsion. For 
radial type of burning appropriate expression for 
I(t) may be used. 

Equation (1) can be written as 

lit) e(f) + [/(f) + K d \ e(0 - A(ty-u(t-t d ) (3) 

from (2), lit) = - [ j 

for the RH-125 roll control system - [K b + /(f)] is 
found to be very small and hence has been neglected 
in the simulation. 

The error function is given by 
e = Q M + m 
The control logic is given by 
u(t) = M for e > 8 

= —M for e <—$ 

= 0 for \e\ <8 t > td 

= 0 for t<t d 

where S — dead zone on the error function. 

M — control torque level. 

The aim of this study is to arrive at ‘optimum’ 
values of the set of parameters K x , K 2 and M. By 
‘optimum’ is meant the maintenance of 0 within 
certain limits demanded by the accuracy require¬ 
ments with less fuel consumption. The simulation 
methods are presented here to achieve the same. 

The best design, in this case would be a system 
with a rate feedback corresponding to the critical 
damping and control torque level and dead-zone of 
such magnitudes that the control torque acts in only 
one direction (opposing the aerodynamic torque) 
throughout the flight time. This avoids unnecessary 
fuel expenditure. 



ANALOG SIMULATION 

The design of roll control system satisfying the 
above criteria is carried out by simulation methods 
which are illustrated here for Rohini-125 rocket. 
The analog simulation discussed here incorporates 
time-varying moment of inertia and the aerodynamic 
roll torque in Fig. 2 3 (approximated as a linear 
function) but it does not include the control delay 
time. The controlled flight time is taken to be 3.6 sec. 

Amplitude Scaling 

Following maximum values are used for amplitude 
scaling of the various quantities : 

8m = 0.10 rad. 

0m = 0.16 rad/sec. 

M m = 0.1 kgm. 

Range of K x and K 2 = 0 to 1.0. 

The values of 0 m and 0™ have been estimated 
from flight test and practical considerations. The 
maximum error corresponding to ^ — 1.0 K 2 — 1.0 
is E m = 0.26 rad. 

Time Scaling 

Time scaling has been done so as to slow down the 
problem by a factor 5 i.e. 

T — 5t where T —Computer time 
t — Problem time 

Results of Analog Simulation 

The simulation has been used to study the effect 
of the following control parameters on the roll angle : 

(1) Control law gains K x , K 2 

(2) Control torque level M 

The following observations are in order : 

(1) Effect of changing for given values of K 2 
and $. Let K t be changed by a factor n. The equa- 
t ; ons of the dead-zone boundaries will then become 
nK$ -j~ ^20 — i S 

Dividing both sides by n, 

KJ + -M = ± B/n. 

n 



Control torque = 0.6 kgm 

K x =1.0; (I)K 2 =0.01, (2) K 2 =0.05, (3) K 2 =0.1, (4) K 2 =0.5 
Fig. 3A. Analog simulation results 
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Fig. 2. Aerodynamic roll inducing torque 
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M=0.1 kgm 

Ki = 1.0, K’2=0.01, K*=0.05, K 2 =0.1, K 2 »0.5 
Fig. 3 B. Analog simulation results 


Hence, the effect of changing K x by a factor n for 
a fixed value of K 2 and 8 will be the same as the 
effect of changing K 2 and S by a factor 1 jn for the 
same value of K 

If both the K x and Ki are changed by a factor n 9 
the effect will be same as the effect of keeping the 

and K 2 unchanged but changing 8 by a factor 1/n. 

(2) The roll angle vs time curves are given in Fig. 3 
for different values of control parameters. The 
curves in Fig. 3 show the effect of various damping 
constant K = (KJKj) on the system for different 
torque levels. The value of K = 0.05 damps the 
system critically and the roll angle 0 is within the 
dead-zone in the smallest time. 

(3) The control torque level of M = 0.5 M m 
where M m is the max. torque level of 0 A kgm, 
is adequate to control the roll angle throughout the 
burning time of 3.6 sec. The aerodynamic roll induc¬ 
ing torque characteristic in Fig. 2 shows that the 


disturbance torque reaches its maximum value of 
0.05 kgm at 3.51 sec., falling sharply thereafter. 
Hence the control torque of 0.05 kgm is adequate 
to control the roll angle throughout the flight. 

Based on the above study the following values of 
control parameters are arrived at : 

8 = 0.05° M = 0.05 kgm. 

Ai = 1.0 K 2 = 0.05 

DIGITAL SIMULATION 

The digital simulation programme incorporates 
time varying moment of inertia, aerodynamic roll 
inducing torque as in Fig. 2 and control actuation 
delay times. Simulation is done for a flight time of 
3.9 sec. 

For the reason given in section 3, variations of 
K 2 , 8 and M are studied for a fixed value of = 1.0. 

DISCUSSION OF RESULTS 

Effect of changing K 2 for given value of K t and 8 

Fig. 4 shows the plot of |0( max. |0| max and 
(T 0 „ x M) observed during the controlled flight 
time of 3.9 secs versus AT 2 . Though the graphs in 
Fig. 4 are given for 8 = 0.4° and M = 0.07 kgm, 
studies were conducted for different values of 8 and 
M and it was observed that the nature of the graph 
was similar. 

It is seen from Fig. 4 that the values of K 2 within 
the range 0.04 to 0.075 are satisfactory from mini¬ 
mum of maximum roll angle deviation, minimum of 
maximum roll rate and minimum control impulse 
point of views. K 2 = 0.05 was selected for studying 
the effect of variation of 8 and M. 
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M -► 

Fig. 5 

Effect of changing 8 and M for given values of K x With these values, the roll history of the system is 
and K 2 shown in Fig. 6. 


After fixing the value of K 2 , various combinations 
of 8 and M were studied, the results of which are not 
given here for brevity. It is observed that for a 
particular value of M the accuracy improves (i.e. 
|0| max decreases) for decreasing value of 8 upto 
a certain value and again deteriorates for further 
decrease in the value of 8. It is also observed that the 
control torque-impulse (T on x M) for a given value 
of M increases continuously for decreasing value 
of 8. By making a trade-olf between accuracy and 
control impulse 8 was selected as 0.2° and from 
Fig. 5, M was fixed as 0.07 kgm. Thus the control 
system parameters are as follows : 

K l = 1.0 K 2 = 0.05 

8 = 0.2° M = 0.07 kgm 



Fig. 6 B. Roll rate vs time 


REMARKS 

Analog simulation was carried out on the AC-20 
Analog Computer (ECIL). The accuracy of 
observations is 0.1% of full scale. Again due to 
the finite operational amplifier gain, the switching 
circuit does not conform to an exact relay type 
nonlinearity, but has a finite slope. The time delay 
in control actuation has been neglected in theanalohg 
simulation. The criterion for parameter variation 
in the analog simulation has been to obtain a critically 
damped system response. In the digital simulation 
(IBM 360/44 system, SSTC), however, the criterion 
for parameter selection is based on fuel optimality 
in addition to the nature of response. As such 
the design values obtained from the digital simula¬ 
tion are more accurate and hence recommended. 

CONCLUSIONS 

The usefulness of simulation techniques in arriv¬ 
ing at the design parameters for a control system have 
been illustrated here. Even though the techniques 
are illustrated for roll control of RH-125, they are 
general in nature and can be applied to any other 
control system. 
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